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Abstract 

We calculate the production of a W boson in association with up to two jets in- 
cluding at least one 6-jet to next-to-leading order (NLO) in QCD at the CERN Large 
Hadron Collider with 7 TeV center-of-mass energy. Both exclusive and inclusive event 
cross section and 6-jet cross sections are presented. The calculation is performed con- 
sistently in the five-flavor-number scheme where both qq' and bq (q ^ b) initiated 
parton level processes are included at NLO QCD. We study the residual theoretical 
uncertainties of the NLO predictions due to the renormalization and factorization scale 
dependence, to the uncertainty from the parton distribution functions, and to the val- 
ues of a s and the bottom-quark mass. 



1 Introduction 

The study of W boson production in association with one and two b jets at both the Fer- 
milab Tevatron collider (pp) and the CERN Large Hadron Collider (LHC, pp) has many 
interesting experimental and theoretical facets. On the experimental side, these processes 
are backgrounds to WH production with the Higgs boson decaying to b quarks, to sin- 
gle top and top-pair production, and to many new physics searches. On the theoretical 
side, these processes offer an interesting testing ground for calculational techniques involv- 
ing heavy quarks with a non-negligible initial-state parton density. Predictions for W boson 
production in association with b quarks are available at higher order in QCD using vari- 
ous calculational techniques (four- flavor [TJ [21 [3J H] and five-flavor number schemes [5]) and 
approximations (massless [HE] and massive b quarks [21 El H]). Recently, NLO fixed-order 
calculations of the qq' — > Wbb parton-level process with massive b quarks have been inter- 
faced with parton-shower Monte Carlo programs within the POWHEG [6] and MC@NLO [3] 
frameworks. 

In this context, the predictions for W + 1 jet and W + 2 jet production with at least 
one b jet include processes where b quarks can have low transverse momentum so that finite 
6-quark mass effects become important. Assuming only massless quarks and gluons in the 
initial state (i.e. working in a four-flavor- number scheme), this signature can only originate 
from the diagram in Fig. [JJ(a), i.e. from qq' — > Wbb, and its higher-order corrections. The 
calculation of NLO QCD corrections to qq' —¥ Wbb of Fig. UJa) with massive b quarks has 
been provided in [2113113], and made available in MCFM [7J. It exhibits interesting theoretical 
features. In particular, large logarithms of the form a s log(mb/p) (where p is a scale of the 
order of the maximum 6-quark transverse momentum) originate from the splitting of a gluon 
into two almost collinear bottom quarks. This happens for the first time in the parton-level 
process qg — > Wb(b)q (where (b) denotes an untagged low p T b quark) depicted in Fig. 121 
This process arises as part of the NLO QCD corrections to qq' — > Wbb, but it is intrinsically 
a tree-level process. As such it exhibits a large renormalization (pn) and factorization (pf) 
scale dependence and, because it is enhanced by large logarithms of the form a s log(mi,/p), 
it potentially introduces a large systematic uncertainty in the calculation, that could be 
tamed only by a complete next-to-next-to-leading (NNLO) calculation of qq' — » Wbb. A 
clever way to reduce this problem is to introduce a 6-quark parton distribution function 
(PDF) [91 [10], defined purely perturbatively as originating from gluon splitting. In this 
way, the scale evolution of the 6-quark PDF resums the large logarithms originating at each 
order and provides a more stable, although approximate solution. In this approach, the 
LO process is considered to be bq — > Wbq', as shown in Fig. [JJ(b), where the 6-quark PDF 
is generated perturbatively from the gluon PDF and the Altarelli-Parisi splitting function 
for g — > bb splitting, and the b is assumed to have too low a pt to be observable. This 
results in exactly the process shown in Fig. [2] with a low-p^ b quark. In this approach, 
the b quark is treated as massless in the hard scattering process qb — > Wbq' , the so-called 
(simplified) ACOT scheme [H [TO] , and its mass only appears as a collinear regulator in 
the initial g — >■ bb splitting function. The resulting logarithms a s Yn.(mb/ Pf) are resummed 
via DGLAP evolution of the 6-quark PDF. The NLO QCD calculation of qb -)• Wbq' has 
been performed in Ref. [5] and made available in MCFM [7]. In fact, as explained in [8], 
the two tree level processes, qq' — > Wbb and qb — > Wbq' and their 0(a s ) corrections can 
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be combined, as long as sufficient care is taken to subtract logarithmic terms that would 
otherwise be double counted. 

In this paper we will combine NLO QCD calculations of qq' —¥ Wbb and qb — > Wbq' 
parton level processes including 6-quark mass effects to provide precise predictions for W + 1 
jet and W + 2 jet production with at least one b jet at the 7 TeV LHC. The choice of the ex- 
perimental signature, jet algorithm, and kinematic cuts has been made according to ATLAS 
specifications [TTJ. We will closely follow Ref. [8] where a consistent combination of these 
two NLO calculations has been performed for the first time to provide predictions for the 
production of a W boson and one 6-jet. It is interesting to note that the calculation of Ref [8] 
has been compared with a measurement of the 6-jet cross section of W boson production 
in association with one and two b jets by the CDF collaboration at the Tevatron [TJ]. This 
comparison found a discrepancy of about two standard deviations [121 H3] • 

After a brief presentation of the theoretical framework in Section (2j we will discuss NLO 
QCD predictions and their residual uncertainties for the 7 TeV LHC in Section [3] and present 
our conclusions in Section HI 




WVV^ 



W 




q' 



W 



(b) 



Figure 1: Leading-order parton- level processes for the production of a W boson and one or 
two jets with at least one b jet. 



2 Theoretical Framework 

The predictions presented in this paper are based on the combination of NLO QCD calcu- 
lations of the qq' — > Wbb [21 |3j |7] and bq -> Wbq [5j parton-level processes, as presented in 
Ref. [8] and implemented in MCFM [7] (where the leptonic W decay is included), and we 
refer to [8] for more details. 

In the NLO QCD calculation of the qq' — > Wbb process the b quark is considered to be 
massive, and only light quarks (q ^ b) are considered in the initial state, i.e. the so-called 
four-flavor number scheme (4FNS) is used. In the NLO QCD calculation of the bq — > Wbq' 
process the 6-quark mass is only kept as regulator of the collinear singularity while it is 
neglected in the hard process so that the hadronic cross section is obtained as follows, 



a b q L ° = J dxidx 2 b(x 1 ,fi) 



^2 Q( x 2, I^F)^b q LO ( m b = 0) + g(x 2 , n F )a^°{m b = 0) 
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W 

Figure 2: A parton-level process contributing to Wbj production that appears at NLO in 
the calculation of 0(a s ) corrections to qq' — > Wbb. This process is also equivalent to the 
LO 6-quark initiated process of Fig. Wljo), with the b quark originating from collinear g —tbb 
splitting. The consistent treatment of this process in the combination of the two NLO 
calculations is described in Section |2j 

An approximate solution of the DGLAP evolution equation for the 6-quark PDF b(x, fip) 
with initial condition b(x, jIf) — at jjlf — m b exhibits the collinear logarithm at leading 
order in a s as follows [HI OB] , 

i {x ,, F) = ^tsl ]0 JtL)r^ Pqa{2}g r±,,\ . (2) 

n \ m bJ Jx z \z J 

When combining the NLO calculation of this process with the NLO calculation of qq' — > Wbb 
this contribution has to be subtracted in order to avoid double counting of the process of 
Fig. [2] which is already included in the 4FNS NLO QCD calculation. The full five-flavor 
number scheme (5FNS) result at NLO QCD, including an all order resummation of collinear 
initial-state logarithms via DGLAP evolution, is then obtained schematically as follows, 

„NLO „NLO f m _l n\ , „NLO 



E 



dx 1 dx2b{x 1 ,n F )q{x 2 ,HF)v b L q O (' m b = 0) ■ (3) 



In fact, the situation is slightly more complicated because the NLO computations of the 
qq' — > Wbb and of the bq — > Wbq' processes are made in two different schemes, one in the 
MS scheme and the other in a decoupling scheme. Hence, in Eq. ([3]) a scheme change is also 
assumed, for which we refer the reader to the literature [9j [TOj IS] for further details. This 
said, we now present the sub-processes relevant for our analysis. In detail, cr^p^ s and cr^ q LO 
in this paper include the following parton level processes: 



qq' — ¥ Wbb at tree level [Fig. [Qa)] and one loop (rrib ^ 0) 
a 4FNS QQ' ~^ Wbbg at tree level (m b ^ 0) 

gq ->■ Wfrfrg' at tree level (Fig. |2]), (m b 7^ 0) 

bq — > Wbq' at tree level [Fig. [Qb)] and one loop (m b = 0) 
a^ q LO bq -)• W^g's 1 at tree level (m b = 0) 

bg —)• Wbq'q at tree level (mj = 0) 

We present results for Vp^f and cr^p^s separately in Tables EH3] for the signatures de- 
scribed in the following section. We leave a full discussion of these results until Section [3j 
but here simply note that the scale dependence of the difference (Jpjff — v^fns clearly shows 
the impact of the initial-state collinear logarithms. The difference is negligible for scales of 
the order of the 6-quark mass but can amount to about 40% of cr^j for /i m 360 GeV. 

Note that we do not include contributions to Wbj production which arise from c — > Wb 
transitions (we assume V cb = 0), since these contributions are suppressed by the smallness of 
V cb and the charm quark PDF. For instance, the dominant contribution to Wbj production 
at the LHC when considering \V cb \ = 0.04 is expected to be eg — > Wbg. Using the setup 
of Table Q] we found for the inclusive W + b event cross section at LO QCD with /1 = /io: 
290 x \V cb \ 2 = 0.5 pb, which is about 1% of the result presented in Table El 

In the calculation of 0"^° in the full 5FNS of Eq. ([3]) we assume the number of light flavors 
to be riif = 5 in both the running of ct s (/iR) and in the determination of the one-loop gluon self 
energy, T> gg , i.e. we only decouple the top quark from the running of a s in the modified MS 
scheme. This choice is motivated by the fact that we usually choose renormalization scales 
fin considerably larger than m?,. Alternatively, one can choose to include the 6-quark mass 
in the calculation of H gg (as done in the implementation of qq' — > Wbb in POWHEG [6]) 
and/or also decouple the b quark (nif = 4). The different treatments generally result in 
differences of about a few per cent in the cross sections presented in Section [3l 

3 Results 

Predictions are provided for W + 1 jet and W + 2 jet production where at least one jet is a 
b— jet as will be measured by ATLAS at the 7 TeV LHC [TT]. Jets are clusters of partons 
built using the anti-fcr algorithm which passed the kinematic cuts specified in Table [U In 
the following, b denotes a jet containing one b quark or one b antiquark, while (bb) denotes 
a jet containing a b and b quark, b and (bb) jets may also contain a light parton. j labels a 
jet without b quarks. We will provide predictions for event and 6-jet cross sections for the 
following signatures: 

• Wb inclusive: one and two-jet events with b jets containing a single b, i.e. Wb + Wbj + 
Wbb. 

• W{bb) inclusive: one and two-jet events with one (bb) jet, i.e. W(bb) + W(bb)j. This 
signature can only result from processes contributing to cr^p^s listed in Section |2j 

• Wb exclusive: one-jet events with one b jet containing a single b, i.e. Wb. 



\Vj\ < 

p= - 


: 2.1 
-l,i? = 0.4 




m b = 
NLO 

a% LC 

Grp = 

v us = 


4.7 GeV 
: CTEQ6.6 [Hj 

\M Z ) = 0.118 
: 1.16639 x 10~ 5 
= V cd = 0.227 


GeV" 2 

(V ub = 



Table 1: Kinematic cuts, jet finding algorithm, PDF sets and input parameters used in this 
study, if not stated otherwise. The kinematic cuts used to simulate the acceptance and 
resolution of the detectors are chosen according to ATLAS specifications [Uj . 

7 TeV LHC: p Tj > 25 GeV 

anti— kx jet algorithm 

M w = 80.41 GeV 

LO: CTEQ6L1 

a L s °(M z ) = 0.130 

gl = 8M W G F /V2 = 0.4266177 

V ud = V cs = 0.974 V us = V cd = 0.227 (V ub = V cb = 0) 



• W(bb) exclusive: one-jet events with one (bb) jet, i.e. W(bb) This signature can only 
result from the processes contributing to a^p^ s listed in Section |2J 

The event and 6-jet cross sections have been obtained consistently at NLO in the 5FNS 
following Ref. [8] and as briefly described in Section [2j 

If not stated otherwise, all results are obtained assuming /xr = /xf = /x, where /xr and \xp 
denote the renormalization and factorization scales respectively. We vary /x between /^o/4 
and 4/xo with \iq = Mw + 2m b . Results for the event cross sections corresponding to the 
four signatures described above are given in Tables [2J131 where we consider non-decaying W 
bosons. The theoretical uncertainty due to the scale dependence can be estimated using these 
results. Inclusive and exclusive event cross sections for pp — > W^bX — > e^vbX assuming 
fi = Ho are provided in Table HJ The results have been obtained by multiplying the total cross 
sections of Tables |2j|3] with the branching ratio BR(iy ± — > e^v) = 0.10805 (labeled as "no 
cuts") and by requiring ATLAS inspired lepton cuts, p e T > 20 GeV, \rj e \ < 2.5, p^ > 25 GeV, 
m^ > 40 GeV, R(l,j) > 0.5 (labeled as 'ATLAS cuts"). All these cuts are implemented in 
the full NLO computation including W decay of MCFM [7] . 

The PDF uncertainties are estimated using the NNPDF2.1 [TF], CTEQ6.6 [16], and 
MSTW08 [19] sets of PDFs as presented in Tables |5]|6j Also shown are predictions for the 
event cross sections for different values of a s obtained with NNPDF2.1 |17J . The dependence 
of our predictions on the value of the 6-quark mass is at the level of a few percent, as can 
be seen from Tables EMI 

The predictions for the event cross sections for W + 1 jet and W + 2 jets with at least 
one b (or (bb)) jet, denoted as cr lj+2 j, W + lb jet (or (bb) jet), denoted as ay, and W + 2 jets 
with at least one b (or (bb)) jet, denoted as a^j-, are provided separately in Table HJ They 
are obtained from the results of Tables EH as follows: 



°ij+2j = [a event (Wb inch) + a cvent (W(bb) inch)} 
<J\j = [a cvent (Wb excl.) + a event (W(bb) excl.)] 

°2j = &1J+2J - CTlj 



Table 2: Inclusive event cross sections (in pb), LHC (-y/i = 7 TeV). No branching ratios or 
tagging efficiencies are included. The Monte Carlo integration error is 0.5%. 





W/+6 incl. 


W + (bb) incl. 


W b incl. 


W-{bb) incl. 




Full 


4FNS 


4FNS 


Full 


4FNS 


4FNS 


fi = /i /4 


66.3 


67.3 


18.6 


40.8 


41.2 


11.4 


/! = /i /2 


60.4 


52.5 


13.8 


37.2 


32.2 


8.6 


^ = ^0 


56.7 


42.6 


10.9 


34.8 


26.3 


6.8 


/! = 2/! 


53.2 


35.5 


8.8 


32.7 


21.9 


5.4 


/i = 4/i 


50.0 


30.1 


7.4 


30.7 


18.7 


4.5 



Table 3: Exclusive event cross sections (in pb), LHC (i/i = 7 TeV). No branching ratios or 
tagging efficiencies are included. The Monte Carlo integration error is within 0.5%. 





W + b excl. 


W+{bb) excl. 


W b excl. 


W~{bb) excl. 




Full 


4FNS 


4FNS 


Full 


4FNS 


4FNS 


fJL = /i /4 


36.7 


36.9 


9.4 


22.8 


22.3 


5.7 


H = ti Q /2 


35.3 


35.2 


7.8 


21.8 


21.5 


4.9 


H = yUo 


33.9 


26.2 


6.7 


20.7 


16.2 


4.3 


fi = 2/xq 


32.2 


22.8 


5.9 


19.8 


14.2 


3.7 


H = ±Hq 


30.3 


19.9 


5.2 


18.8 


12.5 


3.3 



The fc-jet cross sections for W + 1 jet and W + 2 jets with at least one b jet are provided 
separately in Table [TUJ They can be obtained from the Wb and W(bb) inclusive event 
cross sections of Table |2] when the Wbb contribution (normally included in the Wb inclusive 
signatures) is counted twice (since it contains two b jets). More explicitly, using the Wbb 
cross section separately provided in Table [9] in parentheses, the 6-jet cross section can be 
obtained from the event cross sections in Table [2] as follows 

<T fe _jet = [<7 event (W6incl.) - a event (Wbb)] + 2 a event (Wbb) + a event (W(bb) incl.) 

= VeventiWb incl.) + (T event (Wbb) + CT cvent (W(bb) incl.) 

= (?ij+2j + a event (Wbb) 



4 Conclusions 

In this paper we have computed the cross section for the production of a W boson in 
association with up to two jets, including at least one 6-jet, at the 7 TeV LHC. The calculation 
consistently combines next-to-leading order corrections to the parton level processes qq' — > 
Wbb [2J [21 U\ and bq — > Wbq [5j according to the procedure presented in Ref. [S]. We 
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Wb incl. 


W(bb) incl. 


Wb excl. 


W(bb) excl. 




4FNS 


Full 


4FNS 


4FNS 


Full 


4FNS 


W + no cuts 
W + ATLAS cuts 


4.6 
2.2 


6.1 
2.8 


1.2 
0.5 


2.8 
1.3 


3.7 

1.7 


0.7 
0.3 


W~ no cuts 
W~ ATLAS cuts 


2.8 
1.3 


3.8 
1.6 


0.7 
0.3 


1.8 
0.8 


2.2 

1.0 


0.5 
0.2 



Table 4: Inclusive and exclusive event cross sections (in pb), LHC (yfs = 7 TeV), for 
pp — > W^X — > e ± vbX (with fi = fi and CTEQ6.6 [16J). The "no cuts" result is obtained 
by multiplying the total cross sections of Tableland Table[3J respectively, with the branching 

0.10805. The "ATLAS cuts" result is obtained by requiring pf, > 



ratio BR(W /± ->■ e ± u) 

20 GeV, \r] e \ < 2.5, p v T > 25 GeV, 



m 



w 



> 40 GeV, R(l,j) > 0.5. All these cuts are 



implemented in the full NLO computation including W decay of MCFM [7J. 





iy + 6 incl. 


W+{bb) incl. 


W b incl. 


W'(bb) incl. 




4FNS 


Full 


4FNS 


4FNS 


Full 


4FNS 


NNPDF2.1 57J 

ctio ng 

CTEQ6.6 [16] 
MSTW2008 [IS] 


44.1 
42.1 
42.6 
44.2 


59.2 ± 0.7 
56.1 ± 1.5 
56.8 ± 1.4 
59.8 ± 0.7 


11.4 ± 0.1 
10.9 ± 0.3 
10.9 ± 0.2 

11.5 ± 0.1 


27.6 
26.5 
26.3 
28.6 


36.2 ± 0.6 
34.7 ± 1.0 
34.6 ± 1.0 
37.9 ± 0.5 


7.1 ± 0.1 
6.9 ± 0.2 
6.8 ± 0.2 
7.4 ± 0.1 


a s (M z ) = 0.114 
a s (M z ) = 0.124 


39.2 
49.6 


52.0 
66.9 


10.1 

12.8 


24.3 
31.0 


31.6 
41.3 


6.3 

8.0 



Table 5: Inclusive event cross sections (in pb) for different PDF sets including PDF un- 
certainties (NNPDF2.1: full error set. CT10 [HJ/MSTW2008 [T9]: central set), and for 
different values of a s obtained with NNPDF2.1 [17] (with (j,r = (j,f = Ho) ■ 





W + b excl. 


W+{bb) excl. 


W b excl. 


W~(bb) excl. 




4FNS 


Full 


4FNS 


4FNS 


Full 


4FNS 


NNPDF2.1 pi] 

CT10 PS] 

CTEQ6.6 PS] 

MSTW2008 PS] 


26.5 
26.1 
26.2 

27.4 


34.9 ± 0.5 
33.6 ± 1.0 
33.9 ± 0.9 
35.6 ± 0.5 


7.0 ± 0.1 
6.8 ± 0.2 
6.7 ± 0.2 

7.1 ± 0.1 


16.9 
16.3 
16.2 
17.6 


21.6 ± 0.4 

20.8 ± 0.7 
20.8 ± 0.7 
22.4 ± 0.3 


4.4 ± 0.1 

4.2 ± 0.2 

4.3 ± 0.2 
4.6 ± 0.1 


a s (M z ) = 0.114 
a s {M z ) = 0.124 


24.1 
30.4 


31.1 

39.7 


6.2 

7.8 


14.9 
19.0 


18.9 
24.5 


3.8 
4.9 



Table 6: Exclusive event cross sections (in pb) for different PDF sets including PDF un- 
certainties (NNPDF2.1: full error set. CT10 [TB] /MSTW2008 p2]: central set), and for 
different values of a s obtained with NNPDF2.1 [17] (with jjlr = /ip = /xq)- 





W + b incl. 


W + {bb) incl. 


W b incl. 


W~{bb) incl. 




4FNS 


Full 


4FNS 


4FNS 


Full 


4FNS 


m b = 
m b = 


= 4.2 GeV 
= 5.0 GeV 


44.8 
40.8 


58.4 
55.7 


12.8 
10.0 


27.7 
25.2 


35.8 
34.0 


7.9 
6.2 



Table 7: Inclusive event cross sections (in pb) for different values of the b quark mass (m b ) 
obtained with the CTEQ6.6 PDF set [16] and fiR = fip = fiQ. 





W + b excl. 


W+{bb) excl. 


W b excl. 


W~{bb) excl. 




4FNS 


Full 


4FNS 


4FNS 


Full 


4FNS 


m b 
m b 


= 4.2 GeV 
= 5.0 GeV 


27.8 
25.0 


35.3 
33.3 


8.0 
6.1 


17.2 
15.5 


21.5 
20.2 


5.0 

3.8 



Table 8: Exclusive event cross sections (in pb) for different values of the b quark mass (m&) 
obtained with the CTEQ6.6 PDF set [16] and /j,r = [/,f = Mo- 
have particularly focussed on performing a systematic study of our prediction, considering 
a number of sources of theoretical uncertainty, under a set of cuts that will be used by the 
ATLAS collaboration |llj . 

Our results can be summarized as follows, where we have put together all sources of 
uncertainty considered in this paper. The event cross section for W production with one or 
two jets with at least one b jet at NLO QCD at the LHC (7 TeV) as will be measured by 
ATLAS PU is, 



a lj+2j (W+ + W- 



109.2 tlH 



(scale) tH (PDF) ±jj;£ (m b ) pb 



The central value corresponds to the CTEQ6.6 PDF set, // = // and m b = 4.7 GeV. In 
the assessment of the theoretical uncertainty we have considered a very conservative scale 
variation from /xo/4 to 4/xo, i.e. ranging from approximately 20 to 360 GeV. The PDF 
uncertainty is assessed by conmparing the nominal CTEQ6.6 prediction with the results 
obtained for CT10 and MSTW08, while m b is varied from 4.2 to 5 GeV. 
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Table 9: Event cross sections (in pb), LHC (-y/s = 7 TeV) for VT + 1 and W + 2 jet production 
with at least one b jet, W + 1 b jet, and W + 2 jets with at least one b jet (where here b 
jet denotes a jet with a single b or (bb) pair). The event cross sections for Wbb are provided 
separately in parentheses. No branching ratios or tagging efficiencies are included. The 
Monte Carlo integration error is within 0.5%. 





Wu+v 


w+ 


W+ 


Wu + 2i 


w,- 


W 2 ~ 


fi = /i /4 


84.9 [5.6] 


46.1 


38.8 


52.2 [3.2] 


28.5 


23.7 


fi = Ho/2 


74.2 [5.3] 


43.1 


31.1 


45.8 [3.1] 


26.7 


19.1 


^ = ^0 


67.6 [5.0] 


40.6 


27.0 


41.6 [2.9] 


25.0 


16.6 


H = 2/i 


62.0 [4.6] 


38.1 


23.9 


38.1 [2.7] 


23.5 


14.6 


H = 4/^o 


57.4 [4.2] 


35.5 


21.9 


35.2 [2.5] 


22.1 


13.1 



Table 10: 6-jet cross sections (in pb) for W + 1 and W + 2 jet production where at least one 
jet is a b jet, LHC (-y/i = 7 TeV). No branching ratios or tagging efficiencies are included. 
The Monte Carlo integration error is 0.5%. 





°b-jct(W + ) 


(Jb-jct(W') 


fi = /i /4 


90.5 


55.4 


fi = fio/2 


79.5 


48.9 


fi = fi 


72.6 


44.5 


fi = 2fi 


66.6 


40.8 


fi = Afi 


61.6 


37.7 
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